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Abstract. Distant planetary nebulae (PNe) are used to measure distances through the PN lumi- 
nosity function, as kinematic tracers in determining the mass distribution in elliptical galaxies, 
and most recently, for measuring the kinematics of the diffuse stellar population in galaxy clus- 
ters. This article reviews the photometric and spectroscopic survey techniques that have been 
used to detect PNe beyond the Local Group, out to the Coma cluster at 100 Mpc distance. 
Contaminations by other emission sources and ways to overcome them will be discussed as well 
, as some science highlights and future perspectives. 



Keywords, planetary nebulae, surveys, photometry, spectroscopy, distances, elliptical galaxies, 
dark matter, galaxy clusters, intracluster light 



CO 

o 

r-« 
o 

P ■ 1. Introduction 

The search for planetary nebulae (PNe) beyond the Local Group is motivated by 
' several important science projects: 

• The maximum luminosity of a PN in the [OIIL1A5007A line varies by at most 0.1-0.2 
k>( \ mag between the stellar populations of a variety of galaxies. The PN luminosity function 

(PNLF), in particular its bright cutoff, is therefore an important distance estimator (see 
Jacoby 1989, Ciardullo et al. 1989, Jacoby, Ciardullo, Ford 1990, and the review by 
Ciardullo 2003). 

• The radial velocities of PNe, also measured from the redshifted [OIII] line, can be 
used to trace the kinematics of stellar populations that have too low surface brightness for 
other kinematic measurements. PNe are therefore used as tracers of the mass distribution 
in elliptical galaxies at intermediate radii (2 — 5i? e ), where measuring radial velocities 
from integrated absorption line spectra is no longer feasible (Hui et al. 1995, Arnaboldi 
et al. 1996, 1998, Mendez et al. 2001, Romanowsky et al. 2003, Peng, Ford & Freeman 
2004). 

• PNe are currently the only tracers (perhaps excepting globular clusters) that allow 
us to measure the kinematics of the diffuse stellar population in galaxy clusters, and thus 
to investigate the dynamics and cosmological origin of this population (Arnaboldi et al. 
2004, Gerhard et al. 2005). 

PNe have been detected out to distances of some 20 Mpc with narrow-band or slitless 
spectroscopy surveys, and out to 100 Mpc (the Coma cluster) with multislit imaging 
spectroscopy (MSIS) surveys. Sections 2-4 contain short discussions of these survey tech- 
niques. The final Section 5 describes a few recent scientific highlight results based on PN 
surveys and gives a brief outlook on future work. 



1 



2 



Ortwin E. Gerhard 



2. Photometric surveys 

Optical spectra of Galactic PNe show numerous emission lines but essentially no con- 
tinuum flux. The strongest emission line is [OIII]A5007A, which contains ~ 10% of the 
energy output of the central star. Further strong lines are the second [OIII]A4959A line 
(three times weaker than A5007A), and Ha (3 — 5 times weaker). In spectra of PNe at 
Virgo distance, it needs an 8m telescope to clearly see both lines of the [OIII] doublet; 
such spectra only became available recently (Arnaboldi et al. 2003, 2004). 

2.1. Blinking surveys 

Based on the lack of continuum and the strong [OIII]A5007A emission line, Ciardullo, 
Ford, Jacoby and collaborators (e.g., Ciardullo et al. 1989, Jacoby, Ciardullo, Ford 1990) 
developed the so-called on-band/off-band survey technique for distant PNe. This involves 
imaging the field of interest with a narrow-band filter around the [OIII] line and with 
a broad-band filter away from the dominant emission lines. In the early surveys the PN 
candidates were identified by blinking the two images. The PN candidates are those point 
sources in the narrow-band image that have no counterpart in the broad-band image. 
Once identified, the [OIII] fluxes of the PN candidates were determined by aperture 
photometry, and converted to 015007 = — 2.51ogi*5oo7 — 13.74 magnitudes. 

As the early surveys were mostly targeted at elliptical galaxies (in order to derive PNLF 
distances), contamination of the PN samples by HII regions and background galaxies 
was largely absent: the surface density of PNe in elliptical galaxies is much greater than 
the surface density of background galaxies whose redshifted Lya or [Oil] emission lines 
falls into the [OIII] filter. To avoid contamination by Galactic stars, the off-band image 
was taken to be typically ~ 0.25 mag deeper than the on-band image (see also §2.2). 
Nonetheless, at least one such sample contained a substantial fraction of candidates that 
were not confirmed by later work (see Arnaboldi et al. 2003). 

The number density distributions of PN candidates in elliptical galaxies were found 
to approximately follow the galaxy light, indicating constant as, the ratio of PNe to 
B-band luminosity. Inside some radius, of order the effective radius R e , however, the 
number density profile flattens off, as fainter PNe are increasingly lost against the bright 
galaxy background. This effect is present also in more recent slitless surveys in ellipti- 
cals. Construction of the PNLF therefore requires identification of a homogeneous and 
complete PN sample, which excludes the central regions, as well as PNe fainter than a 
magnitude limit set such that the remaining PNe are detected with signal-to-noise of at 
least 9 per pixel (in simulations, essentially all objects are detected at this S/N level). 

2.2. Wide-field surveys 

Following the discovery of three intergalactic PNe in the Virgo cluster near M86 (Arn- 
aboldi et al. 1996), wide-field imaging (WFI) surveys were undertaken to map this much 
more sparse PN population in Fornax, Virgo, and several galaxy groups (Feldmeier, Cia- 
rdullo & Jacoby 1998, Arnaboldi et al. 2002, 2003, Ciardullo et al. 2002, Feldmeier et al. 
2003, Castro- Rodriguez et al. 2003, Aguerri et al. 2005). Again, images are taken in a 
narrow-band filter around the red-shifted [OIII]A5007A line, typically ~ 60A wide to 
cover the wider cluster velocity distribution, and in a nearby broad-band filter. Such 
WFI fields contain many 1000's of stars, so the blinking method is no longer feasible and 
automatic methods for sifting out the PNe must be used. One such technique is based 
on the narrow-band vs. broad-band colour magnitude diagram (CMD), as pioneered by 



Theuns & Warren (1997) An example for the Virgo RCN1 field from Arnaboldi et al. 



(2002) is shown in Figure 1. PNe candidates in this field are those to the right of the 
diagonal line indicating the continuum limit, and below the other lines to exclude Galac- 
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Figure 1. Color-magnitude diagram for all the unresolved sources in the Virgo RCN1 field. 
The horizontal line at m n —m v = — f indicates objects with an observed EW=110 A. The 
diagonal line shows the magnitude corresponding to 1.0 x cr above sky in the V band. Full 
curved lines represent the 99% and 99.9% lines for the distribution of modelled continuum 
objects. The dashed lines represent 84% and 97.5% lines for the distribution of modelled objects 
with m n — m v — — 1. The points are all objects detected by S Extractor. Diamonds ar e objects 



with a significant color excess in the narrow band filter. From Arnaboldi et al. (2002) 



tic stars and background emission galaxies with EW below llOA, as described in the 
figure caption. Because of the much lower surface density of PNe in these intracluster 
fields, the PN catalogues extracted from the images may contain a substantial fraction of 
contaminants. We now briefly discuss these various interlopers and how to avoid them. 

Ordinary Galactic stars. In Fig. 1, stars were excluded from the PN catalogue by 
considering their m n — mi, colour distribution after convolution with the photometric 
errors and discarding all objects above the 99.9% line. However, the photometric errors 
may conspire to result in a bright narrow-band m n and faint broad-band m&, so as to 
result in apparently large EW (negative colour) and undetectable continuum, even for 
stars with AB colour zero (requiring sufficiently large negative error in m n and positive 
error in m?,). Such errors move stars from the top right of Fig. 1 into the region of the 
diagram containing the high-EW sources. While statistically not very likely, for ~ 10 4 
stars in the image this "spill-over" effect (Aguerri et al. 2005) may occur often enough 
to lead to a significant contamination of the PN samples. From the simulations for their 
fields |Aguerri et al. (2005)1 found that 60%/20%/10% of their Virgo intracluster PN 
(ICPN) samples were stars due to the spillover effect, if the limiting magnitude of the 
off-band image was 0/0.5/0.9 mag deeper than the limiting magnitude of the on-band 
image. Thus the contamination by faint stars may be avoided by using sufficiently deep 
off-band images. 

HII regions. Imaging surveys including Ha have found a few compact HII regions 
in one Virgo field (Gerhard et al. 2002, Arnaboldi et al. 2003), probably associated with 
gas that was stripped recently from galaxies, but this is a small fraction of all sources in 
the ICPN catalogue. 
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[Oil] emission galaxies at rcdshift ~ 0.35, whose [Oil] emission lines fall into the 
[OIII] narrow-band filter at Virgo redshift. The contamination from these objects can be 
greatly reduced by retaining only emission sources with observed EW greater than lOOA: 



Hammer et al. (1997) and Hogg et al. (1998) found no [Oil] emitters with restframe EW 
> 70A. This has been confirmed by spectroscopic follow-up observations of Virgo ICPNe 
candidates, none of which showed the [Oil] doublet and a continuum (Arnaboldi et al. 
2004) . Of course, such a stringent EW cut also eliminates some PNe from the catalogue, 
for which the emission is too weak compared to the sky noise. 

Lyman a emitters at redshift ~ 3.1, whose Lya line also comes into the [OIII] 
filter. This class of contaminants was found originally by Kudritzki et al. (2000) taking 
spectra of ICPN candidates in a Virgo field which indeed is consistent with containing no 
intracluster light (ICL; Aguerri et al. 2005, Mihos et al. 2005). Subsequent studies of Lya 
emitters in blank fields have shown that the Lya luminosity function is spatially variable 
and that the brightest Lya emitters are between 0.5 — 1.5 magnitudes fainter than the 
PNLF cutoff in Virgo (Ciardullo et al. 2002b, Castro-Rodriguez et al. 2003, Aguerri 
et al. 2005). With good spatial resolution extended Lya galaxies may be eliminated, 
but according to the simulations of Aguerri et al. (2005) typically ~ 25% of the ICPN 
candidates in Virgo could be Lya background galaxies. 

A large fraction of these will be recognized in follow-up spectroscopy by their contin- 
uum, their broad emission lines, or by the absence of [OIIIJA4959A if this is in the spectral 
range. The spectroscopic results of Arnaboldi et al. (2004) | have confirmed 80% of the 
expected number of PNe in the recent Virgo photometric surveys as real PNe, showing 
that these surveys are now well-understood. Many of the confirmed Virgo ICPNe also 
show the second [OIII]A4959A line besides [OIII]A5007A, confirming without doubt that 
these objects are in the cluster. 

The measured radial velocities then provide us with the kinematics of the intracluster 
stellar population; see the scientific highlight Section [5] Similarly, the combination of 
photometric survey plus spectroscopic follow-up has been used to obtain samples of PN 
radial velocities in the halos of elliptical galaxies. See the references in the Introduction. 



3. Slitless spectroscopy surveys 

3.1. Dispersed/undispersed imaging 

In this technique pioneered by |Mendez et al. (2001) [ undispersed images of the field 
are taken both through a broad band filter and through a narrow band filter centered 
around the redshifted [OIII]A5007A line for the target PNe. The PN candidates are 
identified as those sources that appear only in the on-band image, which also provides 
accurate astrometric positions. Additionally, a dispersed image is taken through a grism 
and the on-band filter. In this image, the PN candidates are still point sources, but are 
now shifted relative to their true positions on the sky by an amount reflecting both the 
mapping by the grism and their radial velocities. By contrast, stars in the field now 
appear as elongated features on the CCD, with length corresponding to the filter width. 

From such data in two overlapping fields, |Mendez et al. (2001) obtained fluxes and 
radial velocities for 535 PN candidates in the E5 galaxy NGC 4697, with velocity errors 
of ~ 40kms _1 and magnitude errors of 0.1 — 0.2 mag for the bright/faint sources. Be- 
cause this is a spectroscopic technique, there are fewer contaminants than in the imaging 
surveys; only background galaxies with faint continuum and narrow emission line will 
remain unrecognized. Mendez et al. obtained a new PNLF distance to NGC 4697 of 
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D ~ 10.5 Mpc, and found that the falling velocity dispersion profile in this galaxy is con- 
sistent with simple isotropic spherical models without dark matter halo out to ~ 2.5R e . 



A similar analysis was done by Teodorescu et al. (2005) for the Fornax elliptical galaxy 



NGC 1344. However, in this galaxy the velocity dispersion profile is not consistent with 
isotropic constant-M/L models, so most likely some dark matter is indicated. For more 
details on this method and the results obtained, see Mendez, this conference. 

3.2. Counter- dispersed imaging (CDI) 

The CDI technique instead uses two images dipersed in opposite directions, to determine 
both the true positions of the PN candidates on the sky and their radial velocities. This 
technique is employed by the special-purpose PN.S instrument mounted at the WHT at 
La Palma, in which the incoming beam, after passing the [OIII]A5007A filter, is split into 
two parts that are simultaneously imaged through separate grisms in two [OIII] arms. 
(A third Ha arm is currently built that will also allow simultaneous Ha imaging of the 
field). Two noteworthy properties of the PN.S instrument are its wide (11. '4 x 10/3) field 
of view and high system efficiency (33%); see Douglas et al. (2002). With the Ha arm, it 
will also be possible to remove high-redshift galaxies with narrow emission line and faint 
continuum, which are the remaining contaminants in the current samples, similar as in 



the Mendez et al. (2001) technique. 



First results based on PN.S data for three intermediate-luminosity elliptical galaxies 



were published by Romanowsky et al. (2003). Like |Mendez et al. (2001)| for NGC 4697, 
they found that all three galaxies (NGC 821, 3379, 4494) have falling dispersion profiles 
out to several R e , which are consistent with dynamical models with little if any dark 
matter out to these radii. The PN.S has also been highly successful in producing a sample 
of some 2500 PNe in M31 (Merrett et al. 2006). Currently, most of the data for the PN.S 
team's core science program have been taken, and are being analyzed with an improved 
data pipeline and revised dynamical models. The goal of the core science program is to 
study a sample of 12 round elliptical galaxies and obtain constraints on their dark matter 
halos out to ~ 6R e (for more details, see Romanowsky, this conference). 



4. Multi-slit imaging spectroscopy surveys 

Beyond about 20 — 30 Mpc distance, PNe are too faint to be detectable with narrow 
band surveys or slitless spectroscopy - their emission disappears in the sky noise in the 
narrow band filter. The brightest PNe in the Coma cluster at 100 Mpc distance have line 
fluxes of 2.2 x 10 -18 ergs -1 cm -2 - this is equivalent to ~ 20 photons per minute through 
the aperture of an 8m telescope, of which ~ 2 will reach the detector for a typical ~ 10% 
overall system efficiency. To detect PNe at such distances requires a spectroscopic blind 
search technique: spectroscopic, so that only the sky noise within a few A dilutes the 
emission from the PN, and blind, because the positions of these faint PNe cannot be 
previously determined. 

The Multi-Slit Imaging Spectroscopy (MSIS) technique (Gerhard et al. 2005) which 
meets these requirements combines a mask of parallel multiple slits with a narrow-band 
filter, centered on the redshifted [OIII]A5007A emission line. Spectra are obtained of all 
PNe that lie behind the slits. The narrow band filter limits the length of the spectra 
on the CCD so that many slits can be simultaneously exposed. For each set of mask 
exposures only a fraction of the field is surveyed; to increase the sky coverage the mask 
can be stepped on the sky. The approach is similar as in some searches for high z > 5 
Lya emitters (e.g., Tran et al. 2004). 

In their pilot study, Gerhard et al. (2005) obtained MSIS observations of the Bernstein 
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Figure 2. Left: Multi-slit mask used in the Coma MSIS observations, superposed on the Coma 
core field. Right: Two-dimensional (top) and one-dimensional (bottom) spectra of 4 PN candi- 
dates and a probable Lya galaxy (right and bottom panel, respectively) in the Coma core field. 
Also shown is a ID arc lamp spectrum at the spectral resolution of th e setup, 7.2A = 440 km s" 1 . 
Note the sec ond [OIII]A4959A line in the third PN spectrum. After |Gerhard et al. (2005)| [See 
online paper "or colour version.' 



et al. (1995) field in the Coma cluster core. This field is at the center of the cluster X-ray 
emission and ~ 5 arcmin away from the cD galaxy NGC 4874. From the Bernstein et al. 
(1995) diffuse light measurement we expected about 425 PNe in this field. The fraction 
of the field surveyed by the mask was ~ 12%; because ~ half of the PNe will be dimmed 
as they are not centered in the slits, we estimated to detect ~ 20 — 30 PNe per mask. 

In the MSIS setup, a compromise must be found between the number of slits and 
hence number of objects detected, and the velocity resolution and signal-to- noise (S/N) 
of the emission sources, which are dependent on the spectral resolution. For the Coma 
observations, a mask with 70 slits was used, each spectrum was 43 pixels on the CCD, 
and the spectral resolution d corresponded to 440 km s -1 . Then the total S/N over all 
pixels becomes 



SNR 



PN 



7.0 x i(r - 4Am x 



-1/2 



1.45A/pix/ V0-"6/ \6<JyJ \3hrs 



-1/2 



1/2 



(4.1) 

for a PN fainter than the cutoff by Am magnitudes. Here <j> is the seeing FWHM, the 
slit width, o y the pixel size, and t exp the exposure time. Relative to the sky noise in one 
pixel, the corresponding S/N of a PN at the PNLF cutoff (Am = 0) was 32. Both values 
assume that the emission from the PN falls through the slit completely. Detection of a 
source is normally considered secure if the S/N relative to the sky noise in one pixel is 
9. In the Coma data this corresponds to SNRpjy ^ 2, i.e., 95% probability of the source 
being real. If this criterion is used for secure detection, PNe about 1.4 magnitudes down 
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the PNLF could be detected. It is clear from equation (4.1) that the most important 
condition for successful detection is good seeing. 

Figure 2 shows the mask used in the Coma observations, as well as the two-dimensional 
and one-dimensional spectra of some of the PN candidates and a probable Lya emitting 
galaxy, and of the arc lamp. The faintest candidate in Fig. 2 has a total SNRp^r = 2.0, 
corresponding to a flux of 3 x 10 -19 ergs -1 cm -2 (~ 68 photons in three hours exposures 
with FOCAS). The reasons to believe that these emission sources are indeed almost all 
PNe in the Coma cluster are: (1) They are unresolved spatially and in wavelength; (2) 
they have no detectable continuum, down to 1.6 x 10~ 20 ergs _1 cm~ 2 A , the flux from 
two O stars at the distance of Coma; (3) for all four bright sources for which [OIII]A5007A 
is sufficiently redshifted, the second [OIII]A4959A line has been detected in the blue part 
of the filter; (4) the emission fluxes are consistent with those expected for PNe in Coma, 
at and below the bright cutoff of the PNLF; (5) the number density and (6) the radial 
velocity distribution are also consistent with the surface brightness of ICL and the radial 
velocities of Coma galaxies in this field. On the contrary, emission sources with measured 
continuum (likely background galaxies) have a much more uniform velocity distribution. 

These first ICPN observations in Coma have already lead the way to the conclusion 
that the two cD galaxies with their associated subclusters are currently undergoing a 
strong interaction prior to merging; see the talk by Arnaboldi in this conference. A lot 
more insight on the dynamics of present-day clusters and their ICL can be expected in 
the next few years - with the MSIS technique, velocities of ICPNe can be measured in 
all clusters out to ~ 100 Mpc distance. 



5. Some recent science highlights and future prospects 

Here I will briefly mention three areas in which recent PN survey work has led to 
important and sometimes surprising results, and where further work will be needed. 

(1) PNe population effects. A reanalysis of the Mendez et al. (2001) sample of 535 
PNe in the elliptical galaxy NGC 4697 has shown that most of the brightest of these PNe 
are part of a subpopulation which is azimuthally unmixed and kinematically peculiar, 
and thus neither traces the distribution of all stars nor can be in dynamical equilibrium 
in the gravitational potential (Sambhus et al. 2006). This introduces an uncertainty in 
the PNLF cutoff luminosity for this galaxy of ~ 0.15mag; the PNLF of two kinematic 
subsamples differ with statistical significance. Probably related is the recent result that 
the bright and faint ICPNe in the MSIS Coma cluster sample (see Sj4| have different 
velocity distributions (Arnaboldi ct al. , in preparation). 

(2) Low density dark matter halos in ellipitical galaxies. The work of Mendez 
et al. (2001) and Romanowsky et al. (2003) based on slitless spectroscopy PN surveys has 
provided evidence that in some intermediate luminosity ellipticals the velocity dispersion 
profiles are falling out to at least several R e , suggesting that the dark matter halos in 
these galaxies must be fairly diffuse. However, other ellipticals appear to have denser 
halos (e.g., Gerhard et al. 2001); also, further study is needed to address the concerns of 



Dekel et al. (2005) The PN.S key project results will be very valuable to clarify these 



issues. 



(3) Kinematics of ICL. Figure 3 based on the results of Arnaboldi et al. (2004) and 



Mihos et al. (2005) gives a vivid illustration how the velocities of ICPNe can shed light 



on the kinematics and dynamics of the ICL, showing clearly the dynamical youth of the 
Virgo cluster. Both M87 and M84 have extended stellar halos in approximate dynamical 
equilibrium with their parent galaxies. Superposed on the M87 halo and in the CORE field 
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Figure 3. Montage of a deep image of the Virgo cluster core from Mihos et al (2005) including 
M 87 (lower left) and M86/ M84 regions (right), and histograms of 1CPN velocities measured 
by Arnaboldi et al. (2004) | in the indicated fields. The image and kinematics show that M87 
has an extended, dynamically bound halo. By contrast, the large halo around M84, M86, and 
NGC 4388 is a superposition; the PNe velocities in this field can clearly be associated with the 
individual galaxies, but M84 also has an extended, dynamically bound halo. The streamers in 
the image as well as the ICPN velocities in the Core and M87 fields show that the intracluster 
stars in Virgo have not yet mixed to virial equilibrium. [See online paper for colour version.] 

stars are seen which clearly have intracluster kinematics but whose velocity distributions 
do not appear to be well-mixed yet, and in the SUB field, all of the PNe observed are 
consistent with belonging dynamically to the large galaxies in the field. ICPN and galaxy 
velocities indicate that, somewhat surprisingly, also the much denser Coma cluster is in a 
stage of strong dynamical evolution (Gerhard et al. 2005, Arnaboldi et al. , in preparation, 
Adami et al. 2005). Further measurements of ICPN velocities as discussed in §4 will show 
whether generally the ICL in nearby galaxy clusters is inhomogeneous and dynamically 
evolving, or whether there are clusters in which it has already found a reasonably mixed 
equilibrium configuration. 

Acknowledgements 

I would like to acknowledge the Swiss National Science Foundation for supporting the 
research program on PNe in distant galaxies and clusters during several grant periods. 



PNe Surveys Beyond the Local Group 



9 



References 

Adami, C, Biviano, A., Durret, F., Mazure, A., 2005, AA, 443, 17 

Aguerri, J.A.L., Gerhard, O.E., Arnaboldi, M., Napolitano, N., et al. 2005, AJ, 129, 2585 
Arnaboldi, M., Aguerri, J.A.L., Napolitano, N.R., Gerhard, O., et al. , 2002, A J, 123, 760 
Arnaboldi, M., Freeman, K.C., Gerhard, O., Matthias, M., et al. 1998, ApJ, 507, 759 
Arnaboldi, M., Freeman, K.C., Mendez, R.H., Capaccioli, M., et al. 1996, ApJ, 472, 145 
Arnaboldi, M., Freeman, K.C., Okamura, S., Yasuda, N., et al. , 2003, AJ, 125, 514 
Arnaboldi, M., Gerhard, O.E., Aguerri, J.A.L., et al. , 2004, ApJ (Letters), 614, L33 
Bernstein, G.M., Nichol, R.C., Tyson, J. A., Ulmer, M.P., Wittman, D., 1995, A J, 110, 1507 
Castro-Rodriguez, N., Aguerri, J.A.L., Arnaboldi, M., Gerhard, O., et al. , 2003, AA, 405, 803 
Ciardullo, R. 2003, in Stellar Candles for the Extragalactic Distance Scale, ed. D. Alloin & W. 

Gieren, Lect Notes Phys., 635, 243 
Ciardullo, R., Feldmeier, J.J., Krelove, K., Jacoby, G.H., Gronwall, C, 2002a, ApJ, 566, 784 
Ciardullo, R., Feldmeier, J.J., Jacoby, G.H., Kuzio de Naray, R., et al. , 2002b, ApJ, 577, 31 
Ciardullo, R., Jacoby, G.H., Ford, H.C. 1989, ApJ, 344, 715 

Dekel, A., Stoehr, F., Mamon, G.A., Cox, T.J., Novak, G.S., et al. , 2005, Nature, 437, 707 

Douglas, N.G., Arnaboldi, M., Freeman, K.C., Kuijken, K., et al. , 2002, PASP, 114, 1234 

Feldmeier, J.J., Ciardullo, R., Jacoby, G.H., 1998, ApJ, 503, 109 

Feldmeier, J. J., Ciardullo, R., Jacoby, G.H., Durrell, P.R., 2003, ApJS, 145, 65 

Gerhard, O., Arnaboldi, M., Freeman, K.C., Okamura, S., 2002, ApJ (Letters), 580, L121 

Gerhard, O., Arnaboldi, M., Freeman, K.C., et al. , 2005, ApJ (Letters), 621, L93 

Gerhard, O., Kronawitter, A., Saglia, R.P., Bender, R., 2001, AJ, 121, 1936 

Hammer, F., Flores, H., Lilly, S.J., Crampton, D., Rola, C, et al. , 1997, ApJ, 481, 49 

Hogg, D.W., Cohen, J.G., Blandford, R., Pahre, M.A., 1998, ApJ, 504, 622 

Hui, X., Ford, H.C, Freeman, K.C., Dopita, M.A., 1995, ApJ, 449, 592 

Jacoby, G.H. 1989, ApJ, 339, 39 

Jacoby, G.H., Ciardullo, R., Ford, H.C, 1990, ApJ, 356, 332 

Kudritzki, R.-P, Mndez, R.H., Feldmeier, J. J., Ciardullo, R., et al. , 2000, ApJ, 536, 19 
Mendez, R.H., Riffeser, A., Kudritzki, R.-P., Matthias, M., et al. 2001, ApJ, 563, 135 
Merrett, H.R., Merrifield, M.R., Douglas, N.G., et al. , MNRAS, in press | |astrc-ph/0603125| ) 
Mihos, J.C, Harding, P., Feldmeier, J., Morrison, H., 2005, ApJ (Letters), 631, L41 
Peng, E., Ford, H.C, Freeman, K.C, 2004, ApJ, 602, 685 

Romanowsky, A. J., Douglas, N.G., Arnaboldi, M., Kuijken, K., et al. , 2003, Science, 301, 1696 
Sambhus, N., Gerhard, O., Mendez, R.H., 2006, A J, 131, 837 

Teodorescu, A.M., Mendez, R.H., Saglia, R.P, Riffeser, A., et al. , 2005, ApJ, 635, 290 
Theuns, T. & Warren, S.J. 1997, MNRAS, 284, Lll 

Tran K.-V., Lilly, S.J., Crampton, D., Brodwin, M. 2004, ApJ (Letters), 612, L89 



